Nature of the Spin Glass State 
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The nature of the spin glass state is investigated by studying changes to the ground state when a 
weak perturbation is applied to the bulk of the system. We consider short range models in three and 
four dimensions and the infinite range Sherrington-Kirkpatrick (SK) and Viana-Bray models. Our 
results for the SK and Viana-Bray models agree with the replica symmetry breaking picture. The 
data for the short range models fit naturally a picture in which there are large scale excitations which 
cost a finite energy but whose surface has a fractal dimension, d s , less than the space dimension d. 
We also discuss a possible crossover to other behavior at larger length scales than the sizes studied. 
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The nature of ordering in spin glasses below the 
transition temperature, T c , remains a controversial is- 
sue. Two theories have been extensively discussed: the 
"droplet theory" proposed by Fisher and Huseld (see also 
Refs. andj-jthje replica symmetry breaking (RSB) 

theory of ParisitJ 0. An important difference between 
these theories concerns the number of large-scale, low 
energy excitations. In the RSB theory, which follows the 
exact solution of the infinite range SK model, there are 
excitations which involve turning over a finite fraction of 
the spins and which cost only a finite energy even in the 
thermodynamic limit. Furthermore, the surface of these 
excitations is arguedH to be space filling, i.e. the fractal 
dimension of their surface, d s , is equal to the space di- 
mension, d. By contrast, in the droplet theory, the lowest 
energy excitation which involves a given spin and which 
has linear spatial extent L typically costs an energy of 
order L e , where 9 is a (positive) exponent. Hence, in the 
thermodynamic limit, excitations which flip a finite frac- 
tion of the spins cost an infinite energy. Also, the surface 
of these excitations is not space filling, i.e. d s < d. 

Recently wdjliil investigated this issue by looking at 
how spin glass ground states in two and three dimensions 
change upon changing the boundary conditions. Extrap- 
olating from the range of sizes studied to the thermo- 
dynamic limit, our results suggest that the low energy 
excitations have d s < d. Siipilar results were found in 
two dimensions by MiddletorO. In this paper, following 
a suggestion by Fisherlij, we apply a perturbation to the 
ground states in the bulk rather than at the surface. The 
motivation for this is two-fold: (i) We can apply the same 
method both to models with short range interactions and 
to infinite range models, like the SK model, and so can 
verify that the method is able to distinguish between the 
RSB picture, which is believed to apply to infinite range 
models, and some other picture which may apply to short 
range models, (ii) It is possible that there are other low 
energy excitations which, are not excited by changing the 
boundary conditionaUta. 

We consider the short-range Ising spin glass in three 
and four dimensions, and, in addition, the SK and Viana- 
Brajc3 models. The latter is infinite range but with a 
finite average coordination number z, and is expected 



to show RSB behavior. All these models have a finite 
transition temperature. 

Our results for the SK and Viana-Bray models show 
clearly the validity of the RSB picture. However, for the 
short range models, our data is consistent with a picture 
suggested by Krzakala and Martir£3 where there are ex- 
tensive excitations with finite energy, i.e. their energy 
varies as L e with 6' = 0, but d s < d. In three dimen- 
sions, this picture is difficult to differentiate from the 
droplet picture where the energy varies as L , because 
of the small value of 9 (~ 0.2, obtained from the magni- 
tude of the change of the ground state energy when the 
boundary conditions are changed from periodic to anti- 
periodicO). It is easier to distinguish the two pictures 
in 4-D, even-though the range of L is less, because 9 is 
much largerta (~ 0.7). 

The Hamiltonian is given by 



n 



(i,3) 



ijSiSj, 



(1) 



where, for the short range case, the sites i lie on a simple 
cubic lattice in dimension d = 3 or 4 with N — L d sites 
(L < 8 in 3-D, L < 5 in 4-D), Si = ±1, and the J l3 
are nearest-neighbor interactions chosen from a Gaus- 
sian distribution with zero mean and standard deviation 
unity. Periodic boundary conditions are applied. For the 
SK model there are interactions between all pairs chosen 
from a Gaussian distribution of width — 1, where 

N < 199. For the Viana-Bray model each spin is con- 
nected with z = 6 spins on average, chosen randomly, 
the width of the Gaussian distribution is unity, and the 
range of sizes is N < 399. To determine the ground state 
we use a hybrid genetic, algorithm introduced by PalE3, 
as discussed elsewheret3. 

Let Sf 5 be the spin configuration in the ground state 

for a given set of bonds. Having found S^°\ we then add 
a perturbation to the Hamiltonian designed to increase 
the energy of the ground state relative to the other states, 
and so possibly induce a change in the ground state. This 
perturbation, which depends upon a positive parameter 
e, changes the interactions Jij by an amount proportional 

toSfsf.i.e. 
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M(e) = e—J2srsPs i S j , (2) 

(*>i> 

where AT;, is the number of bonds in the Hamiltonian. 
The energy of the ground state will thus increase exactly 
by an amount AE^ = e. The energy of any other state, 
a say, will increase by the lesser amount AE^ — e ql°' a ^ , 

where g} '"' 1 is the "link overlap" between the states "0" 
and a, defined by 

^' a) = ^j:S^sfs^S^\ (3) 

in which the sum is over all the Nb pairs where there are 
interactions. Note that the total energy of the states is 
changed by an amount of order unity. 

The decrease in the energy difference between a low 
energy excited state and the ground state is given by 

SE^ = A£<°) - AE^ = e (1 - q^ a) ). (4) 

If this exceeds the original difference in energy, EW - 
E(°\ for at least one of the excited states, then the 
ground state will change due to the perturbation. Wc 
denote the new ground state spin configuration by S^°\ 
and indicate by qi and q, with no indices, the link- and 
spin-overlap between the new and old ground states. 

Next we discuss the expected behavior of q and qi for 
the various models. For the SK model, it is easy to de- 
rive the trivial relation, qi = q 2 (for large N). Since RSB 
theory is expected to be correct, there are some excited 
states which cost a finite energy and which have an over- 
lap q less than unity. According to Eq. (|4|), these have 
a finite probability of becoming the new ground state. 
Hence the average value of q and qi over many sam- 
ples, denoted by [• ■ -] av , should tend to a constant less 
than unity in the thermodynamic limit. This behavior is 
shown in the inset of Fig. |l|. For the Viana-Bray model, 
where there is no trivial connection between q and qi , we 
show in Fig. |l| data for R = (1 — [<fi] av )/ (1 — [?]av) f° r sev- 
eral values of e. This also appears to saturate. We plot 
this ratio rather than [g] av or [qi] av for better compari- 
son with the short range case below. For both models 
we took e to be a multiple of the transition temperature 
(the mean field approximation to it, T^ IF = y/z, for the 
Viana-Bray model) , so that a perturbation of comparable 
magnitude was applied in both cases. 

What do we expect for the short range models? In the 
RSB theory, 1 — [q]av and 1 — [qi] av (and hence the ratio 
R) should saturate to a finite value for large L. To derive 
the prediction of the droplet theory, suppose that the en- 
ergy to create an excitation of linear dimension, I, has a 
characteristic scale of l e (we use 8' rather than 6 to allow 
for the possibility that this exponent is different from the 
one found by changing the boundary conditions). Let us 
assume that large clusters (/ rs L) dominate and ask for 
the probability that a large cluster is excited. The en- 
ergy gained from the perturbation is e(l — qi) ~ e/L( d ~ d ^ 
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FIG. 1. Data for R = (1 - [#]av)/(l - [g]av) for the 
Viana-Bray model with coordination number z = 6 for sev- 
eral values of e. The curvature is a strong indication that 
the data tends to a non-zero value for TV — > oo, as for the 
SK model. The best fits to a + b/N° , shown by the lines, 
give a = 0.872 ± 0.005, 0.883 ± 0.01 and 0.84 ± 0.03 for 
e — \/6, e = y/6/2 and e = \/6/4 respectively. Inset: 1 — [g] av 
for the SK model with the strength of the perturbation given 
by e — 1. Because qi = q 2 , the behavior of [gi] av is very simi- 
lar. The data is clearly tending to a constant at large TV. The 
solid line is the best fit and has a = 0.377 ± 0.004. 

since \/L^ d ~ ds ^ is the fraction of the system containing 
the surface (i.e. the broken bonds) of the cluster. Gener- 
ally this will not be able to overcome the L e energy cost 
to create the cluster. However, there is a distribution of 
cluster energies and if we make the plausible hypothesis 
that this distribution has a finite weight at the origin, 
then the probability that the cluster is excited is propor- 
tional to l / /X d-d »+ 6 ' . In other words 

1 - [<?]av ~ e/i M where fj, = 9' + d- d s . (5) 

As discussed above, 1 — qi is of order 1/ ' L^ d ^ d ^ and so 

1 - [tftW ~ where m = 6' + 2(d - d s ). (6) 

Similar expressions have been derived by Drossel et al.0 
in another context. Eqs. (||) and @ are expected to be 
valid only asymptotically in the limit e — ► 0. In order to 
include data for a range of values of e we note that the 
data is expected to scale as 

1 - fe]av = ^(e/L"), 

l-[qi} av = L-( d - d >)F qi (e/L»), (7) 
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FIG. 2. A scaling plot of the data for [g] av in 3-D according 
to Eq. (0). The data collapse is very good with /i = 0.44. 
The solid curve is a polynomial fit (x 2 = 14.9, d.o.f. = 13), 
constrained to go through the origin, omitting the L = 3 data. 
The dashed line is the linear term in the fit. 



where the scaling functions F q (x) and F qi (x) both vary 
linearly for small x. Note that the above discussion ap- 
plies also to a picture in which 9' — and d s < d. 

A scaling plot of our results for 1 — [g] av in 3D is shown 
in Fig. [2| We consider a range of e from a/6/4 to 4x/6 
(note that T^ IF — a/6) and find that the data collapse 
well onto the form expected in Eq. (Q) with fi — 0.44 ± 
0.02. 

It is also convenient to plot the ratio R, which repre- 
sents the surface to volume ratio of the excited clusters. 
This has a rather weak dependence on e and, as shown 
in Fig. ||, the data for each of the values of e fits well 
the power law behavior L~^ d '~ ds \ expected from Eqs. @ 
and (^]), with d—d s between 0.40 and 0.41 (the goodness 
of fit parameter, Q, is 0.07,0.03,0.85,0.23,0.10, in order 
of increasing e) . The inset to Fig. || shows that there are 
small deviations from the asymptotic behaviour, which 
can be accounted for by a scaling function with the same 
value of fi as in Fig. g and with 

d-d s = 0.42 ±0.02 (3D). (8) 

From this value of /j, and Eqs. (|^) and (||) we find 

9' = 0.02 ±0.03 (3D). (9) 

In order to test the RSB prediction, we tried fits of the 
form R = a + b/L c , which give a = 0.28 ± 0.18, 0.01 ± 
0.14, 0.04 ± 0.11, and -0.28 ± 0.18 (Q = 0.08, 0.01, 0.72, 
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FIG. 3. A plot of R= (1 - te]av)/(l - [g]av), the surface 
to volume ratio of the clusters, in 3D as a function of system 
size for different different values of e. For clarity the e = r/2 
data is omitted. The dependence on e is quite weak, and for 
each value of e the data gives a good fit to a straight line with 
slope [equal to — (d — d s )] consistent with Eq. (^). The inset 
shows a scaling plot of the data according to Eq. ([?]) with the 
same value of fi as in Fig. || The solid curve is a polynomial 
fit (x 2 = 26.0, d.o.f. = 18). 



and 0.52) for e/r = 0.25,0.5,1 and 2. These are con- 
sistent with a = though a fairly small positive value, 
which would imply d s = d, cannot be ruled out. For 
e/r = 4 the fit gives a small positive value, 0.18 ± 0.07 
(Q = 0.79), but this is likely too large a value of e to 
be in the asymptotic regime for these sizes (see the in- 
set of Fig. g). The form R = a + b/L + c/L 2 also fits 
reasonably well the data and gives a between 0.41 and 
0.48 (Q = 0.16,0.03,0.82,0.80,0.16). However, for both 
forms the data are very far from the asymptotic limit 
R ~ a for the sizes considered, unlike for the Viana- 
Bray model (compare the main parts of Figs. ^ and |jj). 
By contrast, the deviation from the asymptotic behavior 
R ~ is quite small (see the inset of Fig. ||). 

In Fig. |i| we show analogous results in 4-D. The calcu- 
lations were performed for two different values e = v8/4 
and a/8 (= T^ IF ). The exponents are essentially the 
same for these two values of the perturbation and the 
fits give fi = 0.26 ± 0.04, d - d s = 0.23 ± 0.02 , and so 
from Eq. (g) we get our main results for 4D: 

9' = 0.03 ±0.05, d - d s = 0.23 ± 0.02 (AD). (10) 

The data in Fig. |] is consistent with the scaling form in 
Eq. (0) but the data for the two values of e are too widely 
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FIG. 4. The main figure shows the ratio 

R — (1 — [<Zi]av)/(l — [g] av ) for different values of e in 4-D. 
The inset shows corresponding data for [q]a,v 



separated to demonstrate scaling. 

Interestingly, our results in both 3-D and 4-D are con- 
sistent with 9' = 0, and, within the error bars, (which are 
purely statistical) incprmatible with the relation 9' = 9, 
since 9 ~ 0.20 in 3-DLM and 9 ~ 0.7 in 4-DEi In 3-D, 
9 — 9' is small, but in 4-D this difference is larger and 
hence the conclusion that 9' ^ 9 is stronger. However, 
the conclusion that d — d s > is less strong in 4-D be- 
cause our value for d — d s is quite small and the range of 
sizes is smaller than in 3-D. 

It would be interesting, in future work, to study the na- 
ture of these excitations to see how they differ from the 
excitation of energy lA-Lsilh 9 > 0) induced by bound- 
ary condition changesElEZnifl. In particular, if their vol- 
ume is space filling, one would expect a non-trivial order 
parameter distribution, P(q), at finite temperatures. 

To conclude, an interpretation of our results for short 
range models which is natural, in that it fits the data with 
a minimum number of parameters and with small correc- 
tions to scaling, is that there are large-scale low energy 
excitations which cost a finite energy, and whose surface 
has fractal dimension less than d. This picture differs 
from the one suggested by Houdayer and MartinE3, in 
which d s = d. Furthermore, the results for short range 
models appear quite different from those of the mean- 
field like Viana-Bray model for samples with a similar co- 
ordination number and a similar number of spins. Other 
scenarios, such as the droplet theory (with 9' = 9 (> 0)) 
or an RSB picture (where 9' = 0,d — d s — 0), require 
larger corrections to scaling, but we cannot rule out the 



possibility of crossover to one of these behaviors at larger 
sizes. 
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